Quasi-two-dimensional bound states 



Jesper Levinsen 1 and Meera M. Parish 1 ' 2 

1 T.C.M. Group, Cavendish Laboratory, J J Thomson Avenue, Cambridge, CBS OHE, United Kingdom 
2 London Centre for Nanotechnology, Gordon Street, London, WC1H OAH, United Kingdom 

(Dated: July 3, 2012) 

We consider the problem of N identical fermions of mass and one distinguishable particle 
of mass mi interacting via short-range interactions in a confined quasi-two-dimensional (quasi-2D) 
geometry. For N — 2 and mass ratios m^/mi < 13.6, we find non-Efimov trimers that smoothly 
evolve from 2D to 3D. In the limit of strong 2D confinement, we show that the energy of the N + 1 
system can be approximated by an effective two-channel model. We use this approximation to solve 
the 3 + 1 problem and we find that a bound tetramer can exist for mass ratios / as low as 5 for 
strong confinement, thus providing the first example of a universal, non-Efimov tetramer involving 
three identical fermions. 



An understanding of the few-body problem can be im- 
portant for gaining insight into the many-body system. 
In dimensions higher than one, few-body bound states 
can, for instance, impact the statistics of the many-body 
quasiparticle excitations. Indeed, for fermionic systems, 
the two-body bound state is fundamental to the under- 
standing of the BCS-BEC crossover [1-4], while the exis- 
tence of three-body bound states of fermions [5, 6] with 
unequal masses can lead to dressed trimer quasiparti- 
cles in the highly polarized Fermi gas [7]. Even in one 
dimension (ID), few-body bound states can impact the 
many-body phase: It has already been shown that one 
can have a Luttinger liquid of trimers [8]. 

In general, attractively interacting bosons readily form 
bound clusters, with the celebrated example being the 
Efimov effect in 3D [9] , where there is a universal hierar- 
chy of trimer states for resonant short-range interactions. 
However, for identical fermions it is less clear whether 
bound states exist, since these are constrained to have 
odd angular momentum owing to Pauli exclusion. For 
short-range interactions in 3D, trimers consisting of two 
identical fermions with mass and one distinguishable 
particle with mass m; can only exist above the critical 
mass ratio m^/m± ~ 8.2 [5], while Efimov trimers only 
appear once m^/m± > 13.6 [10]. However, the existence 
of larger (7V+l)-body bound states involving N > 2 iden- 
tical fermions remains largely unknown — it has only re- 
cently been shown that Efimov tetramers exist in 3D [11]. 

In this Letter, we investigate the problem of N iden- 
tical fermions interacting with one distinguishable par- 
ticle in a confined quasi-2D geometry, where the cen- 
trifugal barrier is reduced and the binding of fermions 
should be favored. Such 2D geometries have recently 
been realised in ultracold atomic Fermi gases [12-16], 
where the fermions are confined to 2D with an effec- 
tive harmonic potential. In addition to allowing one to 
explore the 2D-3D crossover, the harmonic confinement 
can strongly modify the scattering properties of atoms 
via confinement-induced resonances [6, 17, 18]. It has al- 
ready been demonstrated that stable non-Efimov trimers 
can exist for lower mass ratios m^/m± in quasi-2D [6, 19]. 
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FIG. 1: (Color online) Critical mass ratio for the appearance 
of trimers and tetramers in quasi- 2D, where the 2D limit cor- 
responds to tb/uJz — > 0. The solid line follows from the solu- 
tion of the full three-body quasi-2D problem, Eq. (8). Dashed 
lines follow from an effective two-channel model. The vertical 
dotted line marks unitarity, where the 3D scattering length 
diverges. 



Here we show that tetramers involving N = 3 identical 
fermions can appear for m^/m± as low as 5 in quasi-2D 
(see Fig. 1), thus putting it within reach of current cold- 
atom experiments. 

We construct the general equations for the bound state 
of the N + 1 system in quasi-2D and we reveal how to 
simplify the problem in the case of the trimer (TV = 2). 
In the limit of strong 2D confinement, we show that the 
N + 1 problem can be described by an effective two- 
channel model, analogous to that used for Feshbach res- 
onances. This important simplification allows us to solve 
the aforementioned N = 3 problem in quasi-2D. 

In the following, we assume the two atomic species 
{t, 1} to be confined to a quasi-2D geometry by an ap- 
proximately harmonic potential along the z direction, 
Vf^(z) = \m^^J 1 z z 2 . Here, we restrict ourselves to equal 
confinement frequencies for the two species since it allows 
a separation of the relative and center of mass motion 
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along the z-direction, as we discuss below. Such a sce- 
nario can, in principle, be engineered experimentally us- 
ing spin-dependent optical lattices. However, even in the 
case where the confinement frequency is species depen- 
dent, regimes exist in which the few-body properties are 
only weakly affected by this dependence. For instance, 
for large mass ratios and on the molecular side of the Fes- 
hbach resonance, once the t4- dimer is smaller than the 
light atom oscillator length, l\ = ^h/m^uj z , the light 
atom is essentially confined by its interaction with the 
heavy atoms [6]. 

The starting point of our analysis is the T-matrix de- 
scribing the repeated two-body interspecies interaction. 
In the ultracold gases, the interaction is described by a 
zero-range model as the van der Waals range of the in- 
teratomic potential is much smaller than all other length 
scales in the problem, including the confinement lengths. 
The T-matrix may be considered in the basis of the indi- 
vidual motion of a spin-| and t atom. However, due to 
the restriction to equal confinement frequencies for the 
two species, the center of mass and relative motion sep- 
arates and it is advantageous to work in this basis. In 
the center of mass frame of the harmonic oscillator po- 
tential, at energy e below the two-body threshold uj z (we 
set h = 1) and at total 2D momentum q, the T-matrix 
takes the form [20] 



T(q,e) = 



2tt [ I 



-e + q 2 /2(m T + m\) 



where the zero-range interaction is renormalized by the 
use of the 3D scattering length, a s . Here, m r = 
^t m l/( m t + m \) i s the reduced mass and l r z = 
^l/2m r uj z is the confinement length corresponding to 
the relative motion. We use the definition of T [21] 
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The two-dimensional scattering always admits a two- 
body bound state of mass M = + m; and binding 
energy e b > satisfying l r z /a s = T(e h l^ z ). 

The T-matrix in the basis of individual motion is re- 
lated to T by the change of basis 



T n ° ni (a e) 

J- n ni VVL? fc y 



J2 C T ni K.mt)C^ 1 K.'»t) 

n n r n' r 

x^n r (0)^n;(0)T(q,e-n^). (3) 



Here, no and ri\ are the quantum numbers labelling the 
eigenstates of the single-particle Hamiltonians T~L^ = 

\m^J 1 z z§ x while n r and n are the quan- 
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turn numbers in the basis of relative, zqi = z$ — z\, 
and center of mass, Zqi = {tu^zq + m^zi)/M, coordi- 
nates. The wavefunction of the relative motion takes 
the value W(0) = (-l) Wr/ VK - l)H/n r !! if n r is 
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FIG. 2: The diagrams which give the binding energy of the 
TV + 1 bound state in quasi-2D. Black dots indicate the initial 
interaction inside /. 



even, and otherwise. The Clebsch-Gordan coefficients 
Cnnl 1 ( m li m t) = ( n o n i \nn r ) were obtained in Ref. [22] 
and vanish unless no + n\ = n + n r . 

We now turn to the question of the existence of bound 
states consisting of TV spin-^ atoms and a single spin-^ 
atom. To this end, we construct the sum of connected di- 
agrams with TV + 1 incoming atoms (Fig. 2). The t atoms 
are considered on-shell with 2D momenta Iq, harmonic 
oscillator quantum numbers n^, and corresponding single- 
particle energies ek^t — I '^m^-^riiUJ z for i = 1, . . . , TV. 
We consider scattering in the centre of mass frame of the 
2D motion and at a total energy E below the TV + 1 
atom threshold (TV + \)oj z /2. Thus, the 4- atom has 2D 
momentum ko = — k*, harmonic oscillator quan- 

tum number no, and energy Eq = E — ^2^ =1 ek^t- The 
sum of diagrams with TV + 1 incoming particles in which 
the ^ atom interacts first with the t atom numbered 1 



is denoted /^°""^ r . Note that there is no dependence 
on ki as the initial interaction depends only on the total 
momentum of the two atoms. 

The occurrence of a bound state corresponds to a sin- 
gularity of / at its binding energy. This singularity re- 
sults from the summation of an infinite number of dia- 
grams and, at the pole, / satisfies the homogeneous inte- 
gral equation illustrated in Fig. 2: The initial interaction 
is described by a T-matrix, and then the spin-^ atom sub- 
sequently interacts with another of the t atoms. Thus, 
the right hand side contains TV — 1 terms and the integral 
equation satisfied by the bound state energy is (setting 
the volume to 1): 
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where k = ko + ki — k' l5 and the minus sign on the r.h.s. 
appears because / is antisymmetric under the exchange 
of incoming fermions. Equation (4) embodies a simple 
and generic formulation for the (TV + l)-body problem 
in quasi-2D, which in principle allows us to capture the 
crossover from 2D to 3D. Indeed, for the case of TV = 2, 
it is a generalization of the Skorniakov-Ter-Martirosian 
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equation for atom-dimer scattering [23], while for N — 1, 
Eq. (4) simply reduces to the condition for the two-body 
binding energy. Finally, we note that Ref. [24] derived 
an expression similar to our Eq. (4) for the 3D N + 1 
problem. 

An important simplification to Eq. (4) becomes possi- 
ble in the limit of strong quasi-2D confinement, uj z ^> e&. 
Here, the function T can be expanded as 



1 



\n(nx/B) + ^Lx + (D(x 2 ) (5) 
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with B w 0.905 [20, 21]. On the other hand, consider the 
denominator on the r.h.s. of Eq. (4) which we shall write 
for simplicity as e — nu z . Here, the typical energy scale 
e ~ e\y since, for bound states, the function / is strongly 
peaked at momenta ~ v"2m r 65, while it quickly decays 
for large momenta. Now, if we expand the denominator 
in powers of €b/uj z (assuming n ^ 0), then the lowest 
order term vanishes when integrated over momentum due 
to the antisymmetry of /k' .. Consequently, the lowest 
non-vanishing contribution from the denominator is of 
order (e^/co^) 2 when the harmonic oscillator index n is 
non-zero. We conclude that to linear order in €b/uj z the 
integral equation for the N + 1 bound state reduces to 
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with the single particle energies e k = e k o, = E — 
^2iLi e kit' anc ^ T obtained from Eq. (1) using the linear 
expansion of J 7 , Eq. (5). As the effects of confinement 
in this limit are contained solely within the linearized 
T-matrix, Eq. (6) may be obtained through a strictly 
2D 2-channel model [25], where the closed channel corre- 
sponds to excited harmonic oscillator modes. Thus, the 
confinement length V z plays the role of an effective range 
in this model, with the 2D limit l z /a s — >• corresponding 
to a single-channel model. This simplification crucially 
depends on the antisymmetry resulting from Fermi statis- 
tics and it thus does not apply to bound clusters involving 
bosons confined to 2D, as considered in Refs. [26, 27]. Fi- 
nally, we note that a similar simplification was recently 
obtained for quasi- ID atom-dimer scattering [28]. 

We now proceed to solve the three-body problem using 
the above methods. First, note that using the zero-range 
condition and removing the center of mass generally al- 
lows one to reduce the number of harmonic oscillator 
quantum numbers by two in Eq. (4) [6]. For the three- 
body problem, this is achieved by changing coordinates 
to the relative motion of the two atoms initially interact- 
ing, zoi = zq — z\, the relative motion of the pair and the 
third atom, z® 1 = (m^o + m t 2: i)/( 7T H + m t) — ^2, and the 
center of mass Zo 12 = (miZo-\-m^zi-\-m^Z2)/(mi-\-2m^). 
Defining the corresponding quantum numbers noi, n^ 1 , 
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FIG. 3: (Color online) Energy of the trimer in quasi-2D for 
mass ratios m^/m± = 3.5, 4, 5, 6.64 (from top to bottom). 
The solid lines correspond to the full calculation, while the 
dashed lines are derived from the effective two-channel model, 
Eq. (6). 



and .ZVoi2> we adopt the new basis: 
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Here, we have dropped 77,01 from the l.h.s. as the resulting 
equation will be independent of this index. Also, the 
center of mass quantum quantum number, A012, has been 
neglected as it only causes a shift in the energy; since we 
consider the lowest lying trimer, A/012 will be set to in 
the following. In the new basis, Eq. (4) for the trimer 
becomes 

Xkf =T(k 2 ,£-e k2t -n°V) 

02 

^ ^ 02 (0)^ 1 (0)(n° 1 n; )1 |n? 2 no 2 )Xk( 

x > T"n9 \ »( 8 ) 

k' n°tw " £k '^ " £k2t " £k 'i +k ^ " + n02)U)z 

where we can see that the r.h.s. indeed does not depend 
on noi- The matrix element in Eq. (8) may be evaluated 
by a series of coordinate transformations: 
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where several sums can be dropped due to the constraints 
on the Clebsch-Gordan coefficients. 

Since the trimer consists of identical fermions, it must 
necessarily have odd angular momentum L in the x-y 
plane of the 2D layer. Thus, the lowest-energy trimer 
has L = 1, and it can be regarded as a p-wave pairing 
of t fermions mediated by their s-wave interactions with 

n 01 n 01 -a, 

the light I particle. In this case, we have x^l = Xkl e > 
where ^2 is the angle of k2 with respect to the x-axis and 
X is a function of hi = | k2 1 - Integrating over <p2 in Eq. (8) 
then leaves an integral equation that only depends on &2 
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FIG. 4: (Color online) Energy of the dimer (solid line), trimer 
(dotted), and quadrumer (dashed) in 2D as a function of 
mass ratio. The trimer binds when m^/mi > 3.33, consis- 
tent with Ref. [19], while the trimer-tetramer transition oc- 
curs at m^/m^ = 5.0. Inset: The difference between trimer 
and quadrumer energies, Es — E4. 

and nf. The same applies for the two-channel model 
Eq. (6) with N = 2, where now there is only a dependence 
on &2. 

We have calculated the trimer binding energy as a func- 
tion of confinement for a range of mass ratios, as depicted 
in Fig. 3. We see that the binding energy decreases as we 
perturb away from 2D and the centrifugal barrier is in- 
creased. Correspondingly, we find that the critical mass 
ratio m^/m^ for trimer binding smoothly evolves from 
the 2D limit of 3.33 [19] towards the 3D result of 8.2 [5], as 
shown in Fig. 1. For the special case of 6 Li- 40 K mixtures, 
where m^/m^ = 6.64, our results agree with Ref. [6]. In 
the limit of strong 2D confinement, we see that the two- 
channel model captures the lowest-order dependence on 
€b/uj z of the trimer energy and critical mass ratio. 

We can exploit the two-channel model (6) to solve the 
more complicated four-body (TV = 3) problem in quasi- 
2D. Once again, the presence of identical fermions re- 
quires us to consider total angular momentum L = 1. 
Thus, we have for the tetramer 

/k 2 k 3 = /(fc 2 ,fc 3 , A(/> 32 )e^ 2 = -/> 3 ,A: 2 ,-A032)e^ 3 , 

where A(/>32 = ^3 — ^2- We note that a similar equation 
for the tetramer energy was obtained for the 3D problem 
in Ref. [11]. 

Beginning with the 2D limit (e^/uOz = 0), we determine 
the energy of the tetramer compared to the trimer and 
dimer energies (see Fig. 4). Following the transition from 
a dimer to a trimer at mass ratio m^/m^ ~ 3.33, we find a 
trimer-tetramer transition at m^/m^ ~ 5.0. In principle, 
we can use Eq. (4) to consider bound states of even larger 
N, but the problem quickly becomes intractable numeri- 
cally for N > 3. However, we conjecture that composite 
bound states of larger N become possible as m^/m^ is 
increased, since the relative importance of the centrifugal 



barrier between heavy particles (which goes as l/ra^) di- 
minishes compared with the effective attractive potential 
induced by the light particle (~ l/m^j. 

Perturbing away from the 2D limit, we find that the 
trimer-tetramer transition shifts to larger m^/m^ with 
increasing e^/co^, as shown in Fig. 1. Eventually, we ex- 
pect to encounter the four-body Efimov effect in 3D for 
m^/m^ > 13.4 [11]. However, it remains an open ques- 
tion whether our quasi-2D tetramers exist in 3D below 
the critical mass ratio for Efimov physics. 

To conclude, we have provided the first example of 
a universal, non- Efimov tetramer involving three iden- 
tical fermions. Since this quasi-2D tetramer exists for 
mass ratios ra^/ra^ as low as 5, it could potentially be 
probed with ultracold 6 Li- 40 K mixtures. Its small bind- 
ing energy (Fig. 4) suggests that it could appear as a 
resonance in atom-trimer interactions. For instance, if 
a cloud of trimers were prepared under strong quasi-2D 
confinement, signatures of the resonance could be ob- 
served by colliding a cloud of heavy atoms with a cloud 
of trimers, similar to the proposal of Ref. [29] for detect- 
ing an atom-dimer resonance. In addition, the presence 
of trimers and tetramers has implications for the many- 
body phases in quasi-2D, particularly for the highly po- 
larized Fermi gas [30, 31]. 

We emphasize that although we have focussed on the 
N + 1 problem in quasi-2D, the form of Eq. (4) is com- 
pletely general and may be extended to other shapes of 
the confining potential and/or different dimensionalities. 
For instance, in quasi- ID one would use the T-matrix de- 
rived in Refs. [17, 21], along with appropriately redefined 
harmonic oscillator and momentum indices. Futhermore, 
the problem may be studied close to narrow Feshbach res- 
onances, characterized by a large effective range, by using 
an energy-dependent scattering length [32]. Finally, our 
work suggests that a two-channel model may be used to 
model strongly confined quasi-2D Fermi systems in gen- 
eral. 

We gratefully acknowledge fruitful discussions with 
Stefan Baur, Andrea Fischer, Pietro Massignan, Wave 
Ngampruetikorn, and Dmitry Petrov. MMP acknowl- 
edges support from the EPSRC under Grant No. 
EP/H00369X/1. JL acknowledges support from a Marie 
Curie Intra European grant within the 7th European 
Community Framework Programme. 



[1] A. J. Leggett, in Modern Trends in the Theory of Con- 
densed Matter, edited by A. Pekalski and J. Przystawa 
(Springer-Verlag, Berlin, 1980), p. 14. 

[2] C. Comte and P. Nozieres, J. Physique 43, 1069 (1982). 

[3] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 
195 (1985). 

[4] C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, 
Phys. Rev. Lett. 71, 3202 (1993). 



5 



[5] O. I. Kartavtsev and A. V. Malykh, J. Phys. B: At. Mol. 

Opt. Phys. 40, 1429 (2007). 
[6] J. Levinsen, T. G. Tiecke, J. T. M. Walraven, and D. S. 

Petrov, Phys. Rev. Lett. 103, 153202 (2009). 
[7] C. J. M. Mathy, M. M. Parish, and D. A. Huse, Phys. 

Rev. Lett 106, 166404 (2011). 
[8] G. Orso, E. Burovski, and T. Jolicoeur, Phys. Rev. Lett. 

104, 065301 (2010). 

[9] V. N. Efimov, Nucl. Phys. A 210, 157 (1973). 
[10] L. D. Landau and E. M. Lifshitz, Quantum Mechanics 

(Butterworth-Heinemann, Oxford, UK, 1981). 
[11] Y. Castin, C. Mora, and L. Pricoupenko, Phys. Rev. Lett. 

105, 223201 (2010). 

[12] K. Martiyanov, V. Makhalov, and A. Turlapov, Phys. 

Rev. Lett. 105, 030404 (2010). 
[13] B. Frohlich et a/., Phys. Rev. Lett. 106, 105301 (2011). 
[14] A. T. Sommer et a/., Phys. Rev. Lett. 108, 045302 

(2012). 

[15] P. Dyke et a/., Phys. Rev. Lett. 106, 105304 (2011). 
[16] Y. Zhang, W. Ong, I. Arakelyan, and J. E. Thomas, 

Phys. Rev. Lett. 108, 235302 (2012). 
[17] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 
[18] Y. Nishida and S. Tan, Phys. Rev. A 82, 062713 (2010). 
[19] L. Pricoupenko and P. Pedri, Phys. Rev. A 82, 033625 

(2010). 



[20] D. S. Petrov and G. V. Shlyapnikov, Phys. Rev. A 64, 
012706 (2001). 

[21] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
[22] Y. F. Smirnov, Nucl. Phys. 39, 346 (1962). 
[23] G. V. Skorniakov and K. A. Ter-Martirosian, Zh. Eksp. 

Teor. Phys. 31, 775 (1956), [Sov. Phys. JETP 4, 648 

(1957)]. 

[24] R. Minlos, in Proceedings of the Workshop on Singular 
Schro dinger Operators, Trieste, 29 September- 1 October 
1994, edited by G. Dell' Antonio, R. Figari, and A. Teta 
(SISSA, Trieste, 1995). 

[25] E. Timmermans, P. Tommasini, M. Hussein, and A. Ker- 
man, Physics Reports 315, 199 (1999). 

[26] I. V. Brodsky et al, Phys. Rev. A 73, 032724 (2006). 

[27] L. Pricoupenko, Phys. Rev. A 83, 062711 (2011). 

[28] D. S. Petrov, V. Lebedev, and J. T. M. Walraven, Phys. 
Rev. A 85, 062711 (2012). 

[29] J. Levinsen and D. S. Petrov, Eur. Phys. J. D 65, 67 
(2011). 

[30] J. Levinsen and S. K. Baur, arXiv:1202.6564 . 
[31] M. M. Parish and J. Levinsen, in preparation. 
[32] D. S. Petrov, Phys. Rev. Lett. 93, 143201 (2004). 



